The running-mass inflation model and WMAP 
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We consider the observational constraints on the running-mass inflationary model, and in partic- 
ular on the scale-dependence of the spectral index, from the new Cosmic Microwave Background 
(CMB) anisotropy measurements performed by WMAP and from new clustering data from the 
SLOAN survey. We find that the data strongly constraints a significant positive scale-dependence 
of n, and we translate the analysis into bounds on the physical parameters of the infiaton potential. 
Looking deeper into specific types of interaction (gauge and Yukawa) we find that the parameter 
space is significantly constrained by the new data, but that the running mass model remains viable. 



I. INTRODUCTION 

The measurements of the Cosmic Microwave Back- 
ground (CMB) anisotropies provided by the Wilkin- 
son Microwave Anisotropy Probe (WMAP) mission |l| 
have truly marked the beginning of the era of precision 
cosmology. In particular, the shape of the measured 
temperature and polarization angular power spectra 
are in spectacular agreement with the expectations of 
the standard model of structure formation, based on 
primordial adiabatic and nearly scale invariant pertur- 
bations. Assuming this model of structure formation, 
accurate albeit indirect constraints on several cosmo- 
logical parameters have been reported in agreement 
with those previously indicated (see e.g. Q) but with 
much larger error bars. 

Moreover, new, complementary, results from the 
Sloan Digital Sky Survey (SDSS) on galaxy cluster- 
ing (see e.g. |4|)and, more recently, on Lyman- alpha 
Forest clouds (jS]) are now further constraining the 
scenario. 

The question naturally arises if this new cutting- 
edge cosmological data can tell us something about 
inflation, assuming of course that the vacuum fluctu- 
ation of the inflaton generates the primordial pertur- 
bation. (If instead some 'curvaton' field 0, Q does the 
job, the data give information mostly about the prop- 
erties of the curvaton, during some post-inflationary 
era.) It is therefore not a surprise that several recent 
works have attempted to use this new cosmological 
data to constraint and/or falsify inflationary physics. 
In particular, constraints have been placed on the gen- 
eral form of the single field inflationary potential, for 
instance in @, 0, |j3 and @. A different approach 
is to ask about constraints on specific well-motivated 
models, constructed in accordance with present ideas 
about what lies beyond the Standard Model of the 
interactions. A comprehensive survey of such models 
has been provided [TTl IT^ , which in general predict 
an almost scale-invariant spectral index n < 1. Un- 
fortunately, the data do not yet discriminate between 
these scale- invariant models [jf. 

In this paper, we study instead a specific inflation- 



ary scenario, Stewart's running-mass model [T^L Il4t 
ITfSl Hit ITil lisj . a type of inflationary model which 
emerges naturally in the context of supersymmetric 
extensions of the Standard Model. The model is of the 
single-field type (i.e. the slowly-rolling inflaton field 
has only one component [ill Il2|). but nevertheless it 
has the striking signature of a very specific and rela- 
tively strong scale dependence of the spectral index. It 
is therefore a natural question to ask, if such a depen- 
dence is compatible or even preferred by the present 
data: in fact the WMAP collaboration claimed to 
have a slight evidence for non- vanishing running of the 
spectral index in their first year data. Even if their 
best fit value for n' is too strong to be accommodated 
in the usual paradigm of slow-roll models and proba- 
bly is generated mostly by the Lyman a 0] or even 
the low multipole data, we will attempt a conservative 
comparison and repeat our analysis in 20], on order 
to test the predictions of the running-mass model. 

Our paper is organized as follows: In section II we 
discuss the running mass model. In section III we 
present our data analysis method and results. Finally, 
in section IV, we discuss our conclusions. 



II. THE RUNNING MASS MODEL 



A. The inflationary potential 

In common with all supersymmetric models, the 
running mass model [LJ, [lj, lid llfil Il7l Il8| chooses 
for the inflaton cj> a flat direction, in order to sup- 
press all the renormalizable inflaton couplings. Also, 
cj> is many orders of magnitude less than the reduced 
Planck scale Mp = 2.4 x 10 18 GeV which ensures that 
the Planck-suppressed non-renormalizable terms are 
negligible. Symmetries can moreover guarantee that 
odd powers of <f> and in particular the possible lower 
order linear term [2l| are absent. The tree- level po- 
tential is therefore a constant Vq plus a soft super- 
symmetry breaking mass term. Note then that in the 
running-mass model, inflation takes place in a region 
of field space in which one (or more) of the fields are 
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strongly displaced from the vacuum, typically by an 
amount ~ Mp. The model is typically a realization of 
Linde's hybrid inflation [22f. where the displaced field 
is a 'waterfall field', different from the slowly- rolling 
infiaton, and stabilized temporarily during inflation 
by an effective mass term. At tree-level, the running- 
mass model then reduces to the version of hybrid infla- 
tion proposed in [2j| , distinguished from more general 
models by the fact that the waterfall field is displaced 
from the vacuum by an amount of order Mp, and it 
has also a mass not far above the generic supergravity- 
mediated supersymmetry breaking value. 

So in this setting, the infiaton potential is dom- 
inated simply by the soft supersymmetry-breaking 
mass term generated by Vq and its radiative correc- 
tions. These are taken into account by using the 
Renormalization Group (RG) improved potential. To 
construct this, one just needs to substitute the tree 
mass with the running mass 0, Q] : 

V = V + ^m 2 (ln<j))tf + ... . (I) 

Here m 2 (In <f) is obtained by integrating the RG equa- 
tion of the form 

dm 2 dm 2 

dl^ = dhrQ = Pm ' { ) 

with f3 m being the /3-function of the soft infiaton mass 
and depending on all its couplings. A crucial assump- 
tion of the running mass model is that the radiative 
corrections are substantial, as they are exploited to 
realize slow-roll in some region of the potential. In 
fact at the Planck scale <f» ~ Mp, the mass-squared 
is supposed to have the generic supergravity value 
|m 2 | ~ H] = Vb/(3Mp). The radiative corrections 
drive down m 2 , so that, when <f> is many orders of 
magnitude below Mp, it has the much smaller value 
which is needed for viable inflation. There are four 
types of model, depending on the sign of to 2 at the 
Planck scale, and on whether or not that sign has 
changed by the time that the inflationary regime is 
reached. 

Note that sufficiently strong running is realized for 
example in the case of radiative EW symmetry break- 
ing in the Minimal Supersymmetric Standard Model, 
where the running turns one of the Higgs doublet's 
mass from positive to negative: the main difference is 
that in this case it is sufficient to suppress the mass 
and so we do not need to rely on a very large cou- 
pling, as the top Yukawa. Unfortunately, since some 
of the fields take a large Mp v.e.v. after inflation, it 
is not viable to implement the running mass model 
directly within the MSSM, but possibly in some of its 
extensions . 



In general at one loop /3 m is given by 0, ^3 

2C D 
Pm = am 2 + — -dA| 2 m 2 , (3) 

7T l07T z y 

where the first term arises from the gauge interaction 
with coupling a and the second from the Yukawa in- 
teraction A. It is easy to generalize to the case of 
more gauge groups or Yukawas. In the expression 
above, C, D are positive group-theoretic numbers of 
order one, counting the degrees of freedom present in 
the one-loop diagrams contributing to the running, 
m is the gaugino mass, while m 2 oop is the common 
susy breaking mass-squared of the scalar particles in- 
teracting with the infiaton via Yukawa interaction. 
Note that the first term in Eq. © is always nega- 
tive, while the second has no definite sign, since mf oop 
is defined as the mass squared splitting between scalar 
and fermionic superpartners and can have either sign. 
Also the case of a non-interacting infiaton gives di- 
rectly [3 m = and it coincides with the constant mass 
potential. Anyway, in realistic cases the infiaton must 
have some interaction in order to reheat the universe 
or to secure a hybrid end to inflation and so in any 
model we expect naturally some running, even if per- 
haps below the level required by the running-mass 
model 2 . 

Over a sufficiently small range of </>, or for small in- 
fiaton couplings, it is a good approximation to take 
a truncated Taylor expansion of the running mass 
m 2 (ln</>) around a particular scale, which we will 
choose as <fio, the infiaton value at the epoch of hori- 
zon exit for the pivot scale kg; for comparison with 
the WMAP results [|| we choose k = 0.002/i Mpc -1 . 
Then we have: 

V = V Q + im 2 (ln0 o )0 2 - \cH\<? In , (4) 

where we have rescaled the last term w.r.t. 3H] for 
future convenience. The dimensionless constant c is 
proportional to the mass beta function at the partic- 
ularly chosen point, 

/3 m (ln0 o ) 

c = — mT' () 

The linear expansion corresponds on the quantum 
field theory side to the one-loop expansion, since it 
practically neglects the running of B m , which arises at 
two-loops. It has been shown [l|j that for small c, 
as is required by the slow roll conditions, this linear 
approximation is more than sufficient over the range 
of <j) corresponding to horizon exit for astronomically 
interesting scales, i.e. between ko and 8ft. -1 Mpc. 



For a tr ee-level example based on the GUT group SU(6), see 
e.g. 0. 



2 For example one could envisage a very heavy infiaton coupling 
only gravitationally and in that case the running would be 
negligible. Note that we are restricting here to models where 
the /3-function is generated by renormalizable interactions. 
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To simplify the expressions, it is very useful to in- 
troduce a new parameter <f>* via 



m 2 (ln</> ) = 3H]c 



In 



Then Eq. (J7J) takes the simple form 0] 



V = V a --cHjcf> 2 I In | — 



leading to 



y ' * i I 



(6) 



(7) 



(8) 



In typical cases the linear approximation is valid at 
= 0*, and that point is then a maximum or a mini- 
mum of the potential. 

The running-mass model supposes that all soft 
masses at the Planck scale (or some other high scale) 
have magnitude roughly of order 



m 2 , 



L loop\ 



3H 



i • 



(9) 



and that the couplings (gauge or Yukawa) are per- 
turbative, but not too small for the running to be 
substantial. In general then, in the RG equation, 
one product of coupling times mass scale will dom- 
inate over the others and be relevant besides the in- 
flaton mass. Then we expect the coupling c defined 
by Eq. (JSJ to be roughly 



The most economical cases are probably if the two 
quantities are of the same order and/or if the gravitino 
mass during inflation is very near to the vacuum value. 

Taking both assumption seriously, we can get some 
estimate of Hj, depending on the mechanism for su- 
persymmetry breaking: 



Hi 
Hi 
Hi 



10 4 GeV anomaly-mediation 

10 2 GeV gravity-mediation 

10~ 3 GeV gauge-mediation. (12) 



But note that the assumptions above can be easily 
relaxed and any larger value could be plausible. 

In this setting, if the dominant contribution to the 
scalar masses comes from supergravity and barring 
cancellations, we naturally expect Eq. @. For what 
regards the gaugino mass, the expectation of Eq. £J| 
is not so straightforward because gaugino masses can 
be very small with some types of SUSY breaking, de- 
pending on the gauge kinetic function and its depen- 
dence from the moduli fields. But it is not unreason- 
able, as long as the inflationary scale is of the order 
of the gravitino mass in the vacuum. 

We will investigate in the following how good does 
the model compare with the data for reasonable values 
of c and try to reach conclusions on the naturalness 
of the allowed parameters. 



B. The spectrum and the spectral index 



a 



3Hj 



10 



1 to 10 -2 . 



(10) 



|A| 5 



A bigger value of |c| would not allow slow roll inflation 
and infringe upon the validity of our Taylor expan- 
sion Eq. (J2J. On the other hand, using Eq. J7J) as a 
very crude estimate of the inflaton mass at the Planck 
scale, one can see that a much smaller value of |c| is 
probably not viable either, since it would require the 
inflaton mass at that scale to be suppressed below the 
estimate Eq. {jJJ. 

The value of c is directly related to the supersym- 
metry breaking masses in the scalar and gauge sec- 
tor and therefore the theoretical question arises under 
which conditions we expect Eq. © to be satisfied. We 
will assume that supersymmetry breaking originates 
from an F-term not a D-term and that the inflation- 
ary scale Vq is obtained due to the mismatch between 
the two contributions in the SUGRA potential, the 
F-term part and the negative part proportional to the 
superpotential 3|IU| 2 /Mp. Note that we do not need 
to assume that the second term vanishes (so the grav- 
itino remains all the time massive), but it is for our 
purposes sufficient that the F-term during inflation is 
such that 



j j \ infl 

H i > m 3/ 2 ■ 



(11) 



Let us now discuss the more phenomenological is- 
sues of the predicted spectrum and spectral index 
of the primordial density perturbation. We suppose 
that it is generated by the inflaton field perturbation, 
which means that it is purely adiabatic and gaussian. 
It is therefore specified by the curvature perturbation 
1Z(k), with k as usual the comoving wavenumber. This 
quantity is gaussian and hence specified by its spec- 
trum Viz(k)- 

To express such spectrum in the running mass 
model, it is convenient to define yet another parame- 
ter 



cm | 



(13) 



where the subscript denotes as before the epoch of 
horizon exit for the pivot scale fco = 0.002ft, Mpc - . 

Note that s is directly related to a physical param- 
eter of the potential by the relation: 



1 _ m 2 (ln0 o ) 
2 C= "72 



3Hf 



(14) 



so from s and c we can directly access the inflaton 
mass and its couplings at </>q. Note also that the in- 
flaton mass at that scale should be smaller than the 
expected value at Mp and therefore, barring cancel- 
lations, the parameter s is expected to be smaller 
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than unity. We will later plot our results not only 
in the s and c plane, but also in the physical param- 
eters space m 2 4> ) / (iHj) and — /3 m (ln</3 )/(3/fj) 
for fixed value of Hi. 

At the pivot scale, the prediction of the running- 
mass model is 



,1/2 



(*o) 



1 



v,} 



2ti-\/3-Mp|0os 



(15) 

exp(s/c) . (16) 



This quantity is related to the normalization of the 
CMB power spectrum by the relation 8] 



V^ 2 (ka) = 5.43 x 10~ 5 ,4 1/2 = (4.7±0.5) x 10~ 5 



(17) 



where we have used A = 0.751^, obtained by the 
global WMAP fit @ . This normalization can be easily 
satisfied for choices of Vq and <j>o or 0* that correspond 
to reasonable particle physics assumptions. Note that 
cq. I|16l) can also be recast in the form: 



Hi 



27T^ /2 (fc ) 



(18) 



where s also depends logarithmically on cj)Q as given 
in Eq. l(T3|l . This expression can be also used for esti- 
mating |s| 



1 



Hi 



27T^ /2 (fc ) 1 00 | 



(19) 



and therefore the parameter s is also directly related 
to the value of the inflaton field compared to the in- 
flationary scale. So for |s| < 1 we must have 



27T^ /2 (fc ) 



< 



\<fo\ < Mp 



(20) 



where the last inequality stems from the requirement 
of negligible higher order supergravity corrections. 
We see therefore that in this kind of models we ex- 
pect the inflationary scale Hi to be much lower than 
2TrV^ 2 (k )Mp ~ 1.14 x 10 14 GeV. It could be there- 
fore natural to link Vq ~ HfMp to an intermediate 
scale like the supersymmetry breaking scale. Also the 
inflaton value (f>o must indeed be larger than any mass 
splitting and therefore our use of the RGE-improved 
potential is perfectly consistent. 

In this paper we would like to constrain the strong 
scale-dependence of the spectrum, which is given by 



Vn{k) 



— exp 



^ ^cAN(k) _ ^ _ 2cAN{k) 



(21) 

where AN(k) = N - N(k) = ln(fc/fc ). 

To discuss the spectral index and its running, we 
need the first four flatness parameters 0, > given 



by 



1 fMlV"\ 2 s- 

2 i V 



Ml 



^cAN(k) 



M P V " ~ se cAN ( k ) _ c 

M£V'V" 



V 



-cse 



cAN(k) 



,3 — 



M P V V ^ „_2„2cAJV(fc) 

v ~ 



(22) 

(23) 
(24) 
(25) 



The parameters are evaluated at the epoch of horizon 
exit for the scale k. The first parameter e is negligible 
because <fi/Mp is taken to be very small. The con- 
dition for slow-roll inflation is therefore just <C 1, 
which is satisfied in the regime 4> ~ <j>* provided that 
\c\,\s\ <C 1, disregarding the fine-tuned cancellation 
between the two terms. 

Additional and generally stronger constraints on s 
follow from the reasonable assumptions that the mass 
continues to run to the end of slow-roll inflation, and 
that the linear approximation remains roughly valid. 
Discounting the possibility that the end of inflation 
is very fine-tuned, to occur close to the maximum or 
minimum of the potential, we have the lower bound 



> e- cN °\c\ 



(26) 



Note that for negative c, this constraint is very strong, 
requiring a very large value of s even for small c and a 
kind of fine-tuning between s and c to give a reason- 
able value of n — 1 . 

For positive c, we also obtain a significant upper 
bound by setting AiV = Nq in Eq. I|29|l. and remem- 
bering that slow-roll requires | n — 1 1 < 1 : 



s <e 



-cN 



(c> 0) . 



(27) 



In the simplest case, if slow-roll inflation ends when 
n — l actually becomes of order 1, this bound becomes 
an actual estimate, |s| ~ e~ cN °. As discussed in [T^ . 
this upper bound can be relaxed for positive s if the 
running of the mass ceases before the end of slow- 
roll inflation. The approximate region of the s versus 
c plane excluded by these considerations is shown in 
Figure ^ 

Since e is negligible, the spectral index to second 
order is 



i(fc) = l + 277(l + |) +2.13£ 2 



(28) 



This gives 



U{k) ~ 1 = se cAN W (1 _ i.06c) 



l( 



se cAN{k) _ c 



(29) 



The first derivative of the spectral index is given by 



n , (A;) _ dn(fc) ^ 2sce cAN(k)_ 



dln/c 



(30) 
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FIG. 1: The theoretically expected region for the pa- 
rameters s and c for a value of No = 50; the 
solid(red)-line-hatched region is strongly excluded, while 
the dashed(green)-line-hatched region is excluded only if 
the mass is supposed to run up to the end of inflation. The 
dotted (blue) line gives the prediction for the case when 
the end of inflation is triggered by r\ — 1. The circles cor- 
respond to the values in the explicit models discussed in 
section II. C: the upper ones (magenta) refer to the case 
of gauge coupling dominance, while the (blue) one at the 
origin to the case of Yukawa coupling dominance. 



Clearly the spectral index is not constant within cos- 
mological scales unless s or c is very close to zero. 
Note that for this type of models, n' is a higher order 
effect (suppressed by both c and s) as in usual slow- 
roll inflation, but it is not proportional to (n — l) 2 . So 
it is perfectly allowed to have n = 1 at a particular 
scale, with a non-zero running. This is again a con- 
sequence of our initial assumption that in the region 
of the potential where inflation takes place the one- 
loop contribution to the inflaton mass is of the same 
order as the tree one. Note anyway that the perturba- 
tive expansion is not endangered by our assumption, 
since the higher orders stay always smaller than the 
one-loop. Similarly for the slow-roll expansion, the 
second order can become larger than the first one in 
case of strong cancellation between s and c, but the 
perturbative expansion is still solid. 

At the pivot scale, we have, for example, to second 
order in s, c 



no - 1 



2(s- 
2sc. 



c-1.06sc) + -{s-c) 2 (31) 



(32) 



The line s — c/(l — 1.06c) in the s vs c plane cor- 
responds approximately to no = 1. The Harrison- 
Zcldovich case of constant n = 1 is given by the ori- 



gin s = c = 0, while constant spectral index differ- 
ent from 1 is realized either near the c = axis for 
s = (n — 1)/2 or near the s — axis for c = — (n— 1)/2. 

Since the phenomenological parameters only de- 
pend on s — c and sc, as long as higher orders are 
negligible, the allowed region is expected to be sym- 
metric under reflection along the s + c = line. We 
can solve the system of equations exactly and extract 
the parameters s and c from a measurement of no and 
n ; one solution is given as 



1 



ci 



- [no - 1 + 1.06n 



-V( n o-l + l-06n o r-r8n 



1 



si = - [n — 1 + 1.06n 



(33) 



(34) 



-\J{n - 1 + 1.06n ) 2 + 8n 



the second solution is given just by the symmetry, i.e. 
c-2 = —Si and S2 — —c\. 

From the expressions above, it is clear that not all 
values of n' are allowed in the running mass model: 
we obtain the constraint 



n' > 



(no - l) 2 



(35) 



so that a decreasing spectral index is possible only if 
no is different from 1. So in general the prediction of 
the running mass model tends toward positive n', con- 
trary to the result claimed by WMAP 0. Both due 
to this constraint and the exponential dependence on 
c, we see that fitting for arbitrary value of n' is not 
equivalent to performing a fit for the running mass 
model. Note also that, as discussed in [2(|, the fact 
that the running mass model tends in general to give 
more power at low scales for sufficiently large c, can 
give naturally a large value of the reionization redshift 
in the Press-Schechter approximation and accommo- 
date easily the value obtained by WMAP p|. 

We will show in the following the allowed region 
both in the s, c parameter space and in the n' vs no — 1 
plane. 



C. Explicit minimal models 

1. Dominance of the gauge interaction 
The first model to be proposed 0, Q considered 
the case of an inflaton charged under an asymptoti- 
cally free SU(N) group. We will consider the case of 
two matter superfields in the adjoint representation 
as in |17| , where one has very simply a superpotential 
given by 



W = gSTr (<^ 2 ) 



(36) 



then e.g. the direction 
vanishing other fields. 



Note 



(f>\ is D- and F-flat for 
that here for the case 
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of universal SUSY breaking masses, the role of the 
waterfall field can be played by the singlet S and 
the gauge symmetry is unbroken in the true vacuum. 
We can then write easily the /?- functions, 



ft, 



2N 



am 2 



while for the gauge coupling- 
da N 
d\nQ ~~ ~2tt 1 



giving 



a(ln c 



1 



Nc 



■ln< 



(37) 



(38) 



(39) 



So the RG equation for the soft mass can be solved 
analytically to give 



m$(ln0) = ro$(Afp) - 2m 2 (M P ) 
2m 2 (M P ) 
+ (l + fa(Af P )ln# 



(40) 



where we have defined the boundary conditions at 
M P 3 . 

It is then clear that the inflaton mass increases for 
decreasing < Mp and therefore to obtain a reduc- 
tion of the mass absolute value a negative initial mass 
is necessary. 

For the linear approximation, we have then 



2N m 2 (ln0o) 
c = a (In 0o ) — —2 — 

7T Oil j 

27Va(M P ) m 2 (Mp) a 3 (ln0 o ) 
3i?| a 3 (M P ) 

which is a positive number, and 



1 

S + -c 



"^(ln O ) 
3Hj 

mJ(Mp) - 2m 2 (M P ) 

2m 2 (M P ) a 2 (In O ) 
+ 3H] a 2 (M P ) 



(41) 
(42) 

(43) 
(44) 



so s can be positive or negative. 

The power spectrum normalization gives us instead: 



In (0 O ) + ln(27r|s|) = In [V% 2 (k )Hi 



(45) 



The expression above can provide directly an estimate 
of the order of magnitude of 0o, but unfortunately it is 



3 Any function evaluated at Mp is simply the initial value and 
we use ln0 to mean ln(0/Mp). 



not possible to solve directly for this quantity. We can 
instead turn the formula around and use it to define 
the inflationary scale, after we have singled out the 
region where the spectral index is small. 

For 0o we can use as a very rough estimate instead 
the scale where m? exactly vanishes, m= o : 



ln( 



i=0 ) 



-2tt 



Na(Mp) 



1- 1 



K(M P )f 
2m 2 (M P ) 



-1/2' 



(46) 

depending on the value of the coupling constant, the 
scale m =o changes very strongly. Note that in order 
for no to be phenomenologically acceptable, 0o must 
be not much far away: the value of the spectral index 
at m 2 , = is in fact already small, 



n m=Q - 1 



-3c 



(47) 



Assuming a(M P ) ~ 1/24 as in SUSY-GUT models, 
N = 3 and \m 2 AM-p)\ = m 2 (M P ) = 3H], we obtain 
for example 

ro =o = 9.9 x KT 5 Af P ; (48) 

so that the inflationary scale must therefore be 

Hi - 1CT 9 A/ P - 10 9 GeV. (49) 

For the parameters c, s at that point we obtain the 
values: 



c m =o 
s m =o 



0.15 
-0.07 



as expected. This corresponds to 

n{k m=Q ) - 1 = -0.44 



(50) 
(51) 



(52) 



too large to be phenomenologically acceptable, but e.g 
for 0o = 0.25 TO= o one has 



0.16 
0.13 



and therefore 



n(fc ) - 1 = -0-06 
n'(fco) = 0.04. 



(53) 
(54) 



(55) 
(56) 



We see clearly that the spectral index can be very 
small, but with a running of the same order. 

Note anyway that changing a(Afp) by only a factor 
1/2, we obtain much smaller values: 

m =o = 9.74 x 10- 9 A/ P ; (57) 

so that the inflationary scale must therefore be 

Hi - 10~ 13 M P - 10 5 GeV. (58) 
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This strong dependence on the coupling is character- 
istic of dimensional transmutation and allows us to 
construct viable models spanning a very large range 
of Hi, while presenting similar values for the spectral 
index and its running. It is therefore possible to ac- 
commodate practically all the scales in eq. 1121 and in 
general smaller Hi corresponds to smaller gauge cou- 
pling, which for the same gaugino mass means smaller 
c. In fact in this case we have at <pQ = O.250 m= o : 



0.08 
0.06 



and therefore 



n(fco) - 1 = -0.02 
n'(fco) = 0.01 . 



(59) 
(60) 



(61) 
(62) 



2. Dominance of the Yukawa coupling 
We consider here the simplest of the cases studied 
in where the superpotential is given by 



(63) 



and all fields are singlet under gauge interaction. If 
the trilinear susy breaking coupling vanishes, we have: 



dX 
rflnQ 
dm 2 
d\nQ 



3A 
16tt 2 



|A| 2 

3 



(64) 
(65) 



Then for the average scalar mass m 2 = J^. m 2 /3, we 
obtain the simple expression: 



m 2 (ln^/M P ) = 



m 2 (M P 



l-^A 2 (Mp)ln</> 



(66) 



The mass differences instead are constant, so for the 
single masses we have: 

m 2 (ln</>) = m 2 (ln0) - fh 2 (Mp) + m 2 (M P ) . (67) 



m 2 (Mp) + m|(Mp) = 5Hj. Then we have for the 
inflaton mass 



m 2 (ln < 



2H 2 



H 



l-^A 2 (M P )ln0 
and in the linear approximation this gives 
A 2 (0 o )™ 2 (0o) 



(68) 



8tt 2 Hf 



A 2 (M P ) 



1 



12tt 2 

c 2 
2 + 3 



(l-^A 2 (Mp)ln^ ) z 
1 1 
l- 1 f T A 2 (Mp)ln ( /» ~ 2 



(69) 



(70) 



so that c is negative, while s can have either sign. We 
note that in this type of models, the parameter c is 
related to second power of the coupling constant and 
therefore A has to be sufficiently large to give an effect. 

As long as s does not vanish in the interesting re- 
gion, we can find an estimate of <f>Q from the COBE 
normalization. Since the running is slower in this case, 
we can solve iteratively eq. (|43|) for (j) and s, taking 
Hi = 10~ 15 A/ P - 10 3 GeV and A(M P ) = 1. We ob- 
tain then 



M P 



1.3 x 10 



-10 



and 



-0.002 
0.026 . 



(71) 



(72) 
(73) 



We can in this case have |c| < 0.01 since we have 
assumed a somewhat suppressed initial inflaton mass 
m 2 (Mp) = H 2 < 3 iff. So in this case we have very 
small scale dependence: 



n(k ) - 1 = 0.056 
n'(feo) = -0.0001. 



(74) 
(75) 



Note that in this case the masses run from positive to 
negative and that the running is strong at the begin- 
ning and then tends to flatten out at the value given 
by rh 2 = 0. In fact in this case the Yukawa is non- 
asymptotically free and tends to zero at small scales. 
So in the plateau region, if m 2 (Mp) ~ m 2 (Mp), we 
have automatically a flat potential. 

In this specific model, slow roll inflation can be re- 
alized in any of the 3 field directions, but we have 
to consider some non-universal initial masses, since 
the hybrid end is assured only if one of them be- 
comes negative, but not all at the same time. For 
example assume 4>x to be the inflaton and some non- 
universal mass terms coming from higher order term 
in the Kahler potential so that m 2 (Mp) = Hj, but 



III. OBSERVATIONAL CONSTRAINTS 

A. Method 

Our analysis method is based on the computation 
of a likelihood distribution over a grid of precomputed 
theoretical models. 

We restrict our analysis to a flat, adiabatic, A-CDM 
model template computed with a modified version 
of CMBFAST 0, sampling the parameters as fol- 
lows: VLcdmh 2 = u>cdm = 0.05, ...0.20, in steps of 0.01; 
n b h 2 = u b = 0.0018, ...,0.030, in steps of 0.001 and 
h = 0.55, ...,0.85, in steps of 0.05. The value of the 
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cosmological constant A is determined by the flatness 
condition. Our choice of the above parameters is mo- 
tivated by Big Bang Nucleosynthesis bounds on oj^ 
(both from D [25| and 4 He + 7 Li 26] ), from super- 
novae [27} and galaxy clustering observations (see e.g. 

28]). Current data also does not show evidence for ad- 
ditional physics like quintessence (see e.g. [23), extra- 
relativistic particles (see e.g. [13), topological defects 
(see e.g. [31J) or isocurvature perturbations (see e.g. 

32]). We do not consider massive neutrino which may 
have an effect on our results but that are probably 
negligible (see e.g. US]). From the grid above we only 
consider models with age of the universe to > 11 Gyrs. 
We allow for a possible (instantaneous) reionization of 
the intergalactic medium by varying the reionization 
redshift 5 < z r i < 25 and we allow a free re-scaling 
of the fluctuation amplitude by a pre-factor of the or- 
der of C110, in units of C]^q IAP as measured by the 
WMAP satellite. Finally, we let the running param- 
eters c and s vary as follows: —0.2 < c < 0.2, and 
-0.2 < s < 0.2 in steps of 0.008. 

The theoretical models are compared with the 
recent temperature and temperature-polarization 
WMAP data using the publicly available likelihood 
code Q. 

In addition to the CMB data we also consider the 
constraints on the real-space power spectrum of galax- 
ies from the SLOAN galaxy redshift surveys using the 
data and window functions of the analysis of Q . We 
restrict the analysis to a range of scales over which the 
fluctuations are assumed to be in the linear regime 
(k < O^ft-iMpc). When combining with the CMB 
data, we marginalize over a bias b considered as an 
additional free parameter. 

We also include information from the Lyman-alpha 
Forest in the Sloan Digital Sky Survey using the re- 
sults of the analysis of [j| and [3j| which probes the 
amplitude of linear fluctuations at very small scales. 
For this dataset, small-scale power spectra are com- 
puted at high redshifts and compared with the values 
presented in |35j . 



Results 



In Fig. |21 we plot the likelihood contours in the c— s 
plane showing the 1, 2 and 3a contours. The top 
panel is WMAP, middle WMAP+SLOAN and bot- 
tom panel WMAP+SLOAN+Ly-a. As we can see 
there is a strong correlation between the two parame- 
ters along the c — s direction and the inclusion of the 
SLOAN data does not improve significantly the CMB 
constraints. However adding the Lyman-a datasets 
breaks the degeneracy and shrinks the likelihoods. As 
already noticed by |5j , we find that the Lyman a data 
are able to restrict more strongly the scale dependence 
of the spectral index and therefore exclude the param- 
eter space at large |c|; in particular at 68% we have 




w - 



-0.1 - 



0.2 



0.1 



w 



-0.1 



-0.2 



0.2 



0.1 



ui 



-0.1 



-i — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — r~ 



WMAP+SLOAN 




_l I I I I I I I I I I I I I I I I I L 



-0.2 -0.1 0.1 0.2 



l — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — J — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 



WMAP+SLOAN+Ly-a 




_ 2 i— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i— 

-0.2 -0.1 0.1 0.2 



FIG. 2: Likelihood contour plot in the c — s plane 
showing the 1,2 and 3<r contours from the WMAP 
data (Top Panel), WMAP+SLOAN (Middle Panel) and 
WMAP+SLOAN+Ly-a (Bottom Panel). 
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0.08 - 



0.02 - 



0.06 - 



0.04 - 




FIG. 3: Likelihood contour plot in the no — n' plane 
showing the 1,2 and 3u contours from the WMAP 
data (Top Panel), WMAP+SLOAN (Middle Panel) and 
WMAP+SLOAN+Ly-Q (Bottom Panel). 



sc < 0.0043. 

Focusing on the region (s + c > 0), we plot the 
likelihood contours in the 2(s — c) + 1 vs. 2sc plane in 
Figure^ As we explained before, 2(s — c) + l gives the 
value of the spectral index no, while 2sc = n' gives 
the bend in the spectrum. Note that due to the bound 
(|35[) . the viable negative running region is practically 




0.2 



FIG. 4: Likelihood contour plot in the plane c — 
(s + c/2) showing the 1,2 and 3a contours from the 
WMAP+SLOAN+Ly-Q data. We recall that these pa- 
rameters are related to the physical inflaton potential 



parameters by c 
m 2 (In 0o)/(3#f). 



-P m Qn<j>o)/(3Hf) and s + c/2 



indistinguishable from the n' — axis in our scale and 
therefore we do not show it. The maximal negative 
running is in fact —0.0015 at 95% c.l. and it does not 
basically change with the different datasets. As we can 
see from the figure, the WMAP data alone constrains 
n' < 0.05 at 95% C.L.. Including the constraints from 
SLOAN and Ly-a limits the amount of deviation from 
scale invariance to n' < 0.024 at 95% c.l.. 

Finally, in FigQJwe plot the likelihood contours in 
the c vs. s + c/2 plane. As discussed previously, 
this variables correspond to the physical parameters 
in the inflaton potential rescaled by the inflationary 
scale 3Hj. It is clear from the graph that the data 
require a correlation between the inflaton mass and 
the /3-function for large c in order to give a small 
no — 1. It is questionable if such correlation corre- 
sponds to a fine-tuning, and in general depends on 
the explicit realization of the model. In fact in the 
gauge dominated case, we have already emphasized 
that there is relatively large freedom, since the phys- 
ical parameters are more than the observables. For 
the Yukawa dominated case the situation is more con- 
straint, also because the running is weaker and to be 
effective requires always large couplings. In Figure ^ 
are shown the three points in the parameter space that 
we looked at in detail in Section III C. We see that 
WMAP+SLOAN datasets are not able to exclude any 
of the models, but the inclusion of Ly-a excludes the 
first model discussed at 99% c.l.. 
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IV. CONCLUSIONS 

The rather full analysis that we have described con- 
firms the general picture indicated by previous anal- 
ysis [2£j, Ha| . The allowed region in the c vs s plane 
depicted in Figure [21 should be compared with the re- 
gion shown in Figure ^ which approximately delin- 
eates the theoretically disfavoured region, and also 
with the minimum value |c| ~ 10 -2 which is probably 
needed to generate enough running of the mass even 
if we go from the Planck scale to 100 GeV. Combining 
all of these, we see that if |c| is significantly above the 
minimum value, only the version of the model with c 
and s both positive is viable. In that case, the spectral 
index has significant running which will be detectable 
in the foreseeable future. On the other hand, if |c| is 
really of order 10 -2 , all choices of the signs of c and s 
are possible except maybe negative c with positive s. 
Furthermore, if that extreme case can be realized in a 
viable running-mass model the running of n will be so 
small that it may never be detectable. The data are 
now starting to squeeze the allowed region to values 
\c\ < 0.1, still away from the lower bound 10 -2 . The 
smallness of |c| could be interpreted as a hint that the 



inflationary scale needs probably to be low, to make 
the running from Mp effective. Note anyway that 
Hj values of the order of 100 GeV, the expected soft 
SUSY breaking masses in the true vacuum from grav- 
ity mediation, are still acceptable, as demonstrated in 
the cases of the simple models presented. 

Looking at the observational situation in more de- 
tail, our results show again that the CMB data can 
put very strong constraints on the value of the spec- 
tral index at large scales, no = 1 + 2(s — c), but still 
allow a pretty large scale-dependence. Other infor- 
mation on the power spectrum, like Lyman a data, 
are needed to reduce the parameter space in the or- 
thogonal direction. Even with this inclusion, though, 
values of |c| of the order of 0.1 are allowed and we 
have n' < 0.037 at 99% c.L. Our allowed region also 
looks still symmetric under reflection with respect to 
the s + c = line: this means that the present data 
are not sensitive enough to distinguish the variation 
of n' that is predicted by the running-mass model. 
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